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Abstract 


We study quasi-Monte Carlo integration in a weighted anchored Sobolev space. 

As the underlying integration nodes we consider Halton sequences in prime bases 
p = (pi,... ,p s ) which are shifted with a p-adic shift based on p-adic arithmetic. 

The error is studied in the worst-case setting. In a recent paper, Hellekalek together 
with the authors of this article proved optimal error bounds in the root mean square 
sense, where the mean was extended over the uncountable set of all possible p-adic 
shifts. Here we show that candidates for good shifts can in fact be chosen from a 
finite set and can be found by a component-by-component algorithm. 
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MSC: 65D30, 65C05, 11K38, 11K45. 

1 Introduction 

We study the problem of approximating the value of the integral I s (f) Jj 0 1 j s f{x) dec of 
functions / belonging to a reproducing kernel Hilbert space 'H(K) of functions [0, l] s —> R. 
One way of numerically approximating I s (f) is to employ a quasi-Monte Carlo (QMC) 
rule, 



where P/v,s = { x o , *i> • • • , x n-i} is a set of N deterministically chosen points in [0, l) s . 
It is well known (see, e.g., [31SUSUIIUI3]) that point sets which are in some way evenly 
distributed in the unit cube yield a low integration error when applying a QMC rule for 
approximating /,(/). 

We study the error of QMC rules in the worst-case setting. The worst-case error of 
an algorithm Qn, s based on nodes Pn.s is defined as the worst integration error over the 
unit-ball of 'H(K), i.e., 


e N , s { p N,sJ<) = sup \I a (f) - QnM)\- 


fen(K) 


*The authors are supported by the Austrian Science Fund (FWF): Projects F5506-N26 (Kritzer) and 
F5509-N26 (Pillichshammer), respectively, which are part of the Special Research Program ’’Quasi-Monte 
Carlo Methods: Theory and Applications”. 
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An essential question in the theory of QMC methods is how the sample nodes Pn, s of a 
QMC rule Qn,s should be chosen. 


Shifted Halton sequences. In this paper we focus on a special kind of point sequences 
underlying a QMC rule, namely Halton sequences (cf. [5]) whose definition is based on 
the radical inverse function. Let p > 2 be an integer, M = {1,2,3,...}, and N 0 = NU{0}. 
For n £ No, let n = + npp + ri 2 p 2 + ■ ■ ■ be the base p expansion of n (which is of 

course finite) with digits n, £ {0,1,... ,p — 1} for i > 0. The radical inverse function 
4> p : N 0 —> [0,1) in base p is defined by 


M n ) '■= 

r =0 


n r 

pr+l • 


Halton sequences can be defined for any dimension s £ N. Let pi,... ,p s > 2 be s integers, 
and let p = (j)\.....p s ). Then the s-dimensional Halton sequence H p in bases pi,... ,p s 
is defined to be the sequence H p = ( x n ) n > 0 C [0, l) s , where 


x n = (0 Pl (n),0p 2 (n),... ,0 Ps (n)), for n £ N 0 . 


It is well known (see, e.g., ID]) that Halton sequences have good distribution properties 
if and only if the bases pi,...,p s are mutually relatively prime, and for the sake of 
simplicity we assume throughout the rest of the paper that p = (p\ ,... ,p s ) consists of s 
mutually different prime numbers. 

We also introduce a method of randomizing the elements of the Halton sequence which 
is referred to as a p-adic shift. This special case of randomization is based on arithmetic 
over the p-adic numbers and is perfectly suited for Halton sequences H p . 

Let p be a prime number. We define the set of p-adic numbers as the set of formal 
sums 

{ OO 

z = ^2 z rP r '■ z r ^ {0,1, • • •, p — 1} for all r £ N 0 

r=0 

Clearly No C Z p . For two nonnegative integers y, z £ No C Z p , the sum y + z £ Z p is 
defined as the usual sum of integers. The addition can be extended to all p-adic numbers. 
The set Z p with this addition, which we denote by +^ p , then forms an abelian group. 

As an extension of the radical inverse function defined above, we define the so-called 
Monna map 

OO 

(j) p :Z p ^ [0,1) by <j> p (z) := ^ (mod 1) 

r=0 ^ 

whose restriction to No is exactly the radical inverse function in base p. In order to keep 
the used notation at a minimum we denote both, the Monna map and the radical inverse 
function, by (j) p . We also define the inverse 



0+ : [0,1) ->• Z p by 0+ 



5 ^x r p r , 

r=0 


where we always use the finite p-adic representation for p-adic rationals in [0,1). By a 
p-adic rational, we understand a number in [0,1) that can be represented by a finite p-adic 
expansion. 

For a prime number p and for x £ [0,1) we consider the following p-adic shifts: 
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• p-adic shift: for a G [0,1), we define x © p cr G [0,1) to be 

x® p a = 0 p ( 0 j(a:) +z p <t>p{o))- 

• simplified p-adic shift: for m G N and a G [0,1), we write x cr to be the 
truncation of x © p cr to the m most significant digits, i.e., if 4>p(x) +z p 0 p cr = 

VrV r ~ X e Z p , then 


x ©pZ a = fa \Y VrP 


,r—l 


, r=l 


• mid-simplified p-adic shift: for m G N and a G [0,1), we write 


1 

2p m ’ 


If the choice of m is clear from the context, we may often omit m in the notation ffi® 1 ^ 
and ffi™^ and write © p mp and ffi™ ld instead. 

In the s -variate case, for given bases p = (pi,... ,p s ), a point x = (xi,..., x s ) G [0, l) s , 
and given <r = (<ti, . .., cr s ) G [0, l) s and m = (mi ,..., m s ) G the above shifts are 
defined component-wise and we write x(B p cr G [0, l) s , xffi^P cr and x&p^cr, respectively. 

For a point set Y = {y n : n = 0,..., iV — 1} we write 

Y © cr := {y n © cr : n — 0,..., TV — 1} where © is either © p , ffi p n ^, or ffi™ 1 ^ . 


A weighted Sobolev space. In this paper, we are going to consider the problem of 
numerical integration of functions / that belong to a weighted anchored Sobolev space. 
Before we give the definition we introduce some notation which we require for the fol¬ 
lowing: assume that 7 = ( 7 J )°fi 1 is a non-increasing sequence of positive weights, where 
1 > 71 > 72 > ■ ■ •. These weights are used in order to model the influence of the different 
variables of the integrands, an idea which was introduced by Sloan and Wozniakowski 
[T5] . For s G N let [s] := {1,..., s}. For u C [s], x u denotes the projection of x G [0, l] s 
onto [0, l]l u ! consisting of the components whose indices are contained in u. Furthermore 
we write (x u , 1) G [0, l] s for the point where those components of x whose indices are not 
in u are replaced by 1 . 

We consider a weighted anchored Sobolev space T-L(K sn ) with anchor 1 = (1,1,..., 1) 
consisting of functions on [ 0 , l] s whose first mixed partial derivatives are square integrable. 
This space is a reproducing kernel Hilbert space with kernel function 

S 

K s ,j(x,y) = Y[(l + ^ j mm(l-x j ,l-y j )) for x, y G [ 0 , l] s , ( 1 ) 

j =1 

where x = (aq, x 2 , • • •, x s ) and y = (pi, y 2 ,..., y s )- The inner product is given by 

f <9l u l <9l u l 

(/> 9)k.„ = Yj W 1 / /'(aw l) T —g(x u , 1) dx u . 

c[ 1 ^[0,l] |u| OX u OX » 
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Here 7 U = Iljeu 7jj i 11 particular 70 = 1. Furthermore, we denote by §^~h the derivative 
of a function h with respect to the x 3 with j e u. The norm in is given by 

ll/lliq,„ = V (/> f)K an - The Sobolev space 'H(K sr/ ) has been studied frequently in the 
literature (see, among many references, e.g. m 0 ei ei uni ca H5i usd. 


It is well known that the squared worst-case integration error in a reproducing kernel 
Hilbert space can be expressed in terms of the kernel function. In the particular case 
of the kernel K S)1 , it is easily derived with the help of [3] Proposition 2.11] that for 
Pn, s = {^ 0 , • • •, apv-i} bi [0, l) s , where x n = (x nA ,..., x UjS ) for n = 0,1,..., N - 1, we 
have 


'N,s 


(p N ,s,K sn ) = nf 1 


3 


N-l s 


N 


1=1 n =0 i= 1 

_ N-l s 




| n-xi 


xi n i 1+7j ““i 1 - x “. j ' 1 - Xh <^ ■ 


( 2 ) 


n,h =0 i=l 


Hence the worst-case error can be computed at a cost of 0(sN 2 ) arithmetic operations. 


In [ 6 ] the authors studied the root mean square worst-case error in H{K srf ) of the 
p-adic shifted Halton sequence extended over all p-adic shifts, i.e., 

e N ,s(H p , K S1 ) : = y^E a[e 2 NjS {H p ® p rr. K sn )]. 


The following result is the main result of [ 6 ]. 
Theorem 1 ( [ 6 ], Theorem 1]). Let N > 2. We have 


[^n,s(H p , K sn )] < 


-i=i 


I (i + y^loglV) 


Pi 


log p. 


n 

3 = 1 


i + £ 
2 


n 

3 = 1 


i+ 

6 J 


(3) 


In particular, if YffpLi Ij 777 “ < 00 > ^ en f or an y d > 0 we have 

e N,s{Hpt K s ~f) <^5 j7iP ’ 

where the implied constant is independent of the dimension s. 

The bound ([3]) is, up to log-factors, optimal. For a further discussion of the result, 
especially with respect to the dependence on the dimension s we refer to [6]. Theorem Q] 
can also be interpreted in the “deterministic” sense that for every fixed N > 2 there 
exists a p-adic shift er e [0, l) s such that the squared worst-case error of the initial N 
elements of the corresponding p-adically shifted Halton sequence satisfies the bound (J3J). 
The problem with this interpretation is that the p-adic shift has to be chosen from an 
uncountable set, namely the s-dimensional unit cube. This is a big drawback if one wants 
to effectively find good p-adic shifts. 

It is the aim of this short paper to show that it suffices to choose the p-adic shifts, which 
yield an upper bound of the form (j3J), from a finite set. This set of possible candidates 
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has size N s which is of course huge already for moderately large s or N. However we also 


show, that in principle good shifts can be found by a component-by-component (CBC) 
algorithm. This idea is borrowed from the construction of good lattice point sets which 
goes back to Korobov [7] and to Sloan and Reztsov [13], and which is nowadays used in 
a multitude of papers. With this “adaptive search” the search space is only of a size of 
order 0{sN). 

The rest of the paper is structured as follows: In Section [2] we prove some auxiliary 
results. The CBC construction of p-adic shifts as well as the statement and proof of the 
main result of this paper are presented in Section [3] 

2 Auxiliary results 

We use the following notation: for p e N and m G N 0 let 


Q (P m ) ■■= {ap~ m : a = 0,1,... ,p m — 1}. 


We now show the following lemma. 

Lemma 1. Let H Pt N be the point set consisting of the first N elements of H p and let 
m G N be minimal such that N < p m . Furthermore, let a m e Q{p m ). Then it is true that 



Proof. Let H Pj n = {h 0 , hi,, /ijv-i}- From fl2]) we obtain 


l 


'P 



P 


m 




n=0 
AT— 1 



n,k =0 
n^k 


1 + 7! min {1 - (K ® p (a m + 5)), 1 - (h k ® p (a m + h))} d<5. 


For given n G {0,1,..., — 1}, let us now analyze the quantity 


hn ©p (©7i + *5) {fin) +Z P fi p {&m + &)) • 


(r) 

The base p expansion of h n is of the form h n = ^+, since N < p m . Furthermore, 

the base p expansions of a m and 5, respectively, are of the form 


oo x(r) 

and 




5 


r=m-\-1 


due to the assumptions on a m and 6. Consequently, 

m m oo 

<f>}(hn) = X^SrV -1 and = 5>hy-y Zp <j (r y _i . 

r= 1 r=l 

Let 

m+1 

<f>i(hn) +Z P <f>p{c m ) = Y1 VrP’'- 1 

r=1 

with y r G {0,1,. .. ,p — 1}. Then we obtain 

m 

4>p ( K ) + Zp (Pp {c m + S) = Y Vrf~ X +Z P ym+lP m +Z P <f>f (5) , 


r=1 


Note that Y1T=i VrP r 1 is the truncation of the p-adic sum 0+ {h n ) +z p 4>p{c m ) to the first 
m digits. Hence 


Y PrP r 1 ) = hr 


Nsmp 




, r=1 


For short we write 


^(/i n ,(T m ) . 4 ) p{lJm+lP )• 

Note that 0+ (£(h n , cr m )) = y m+ ip rn . Hence we can write 

K ® P {c m + 8) = <Pp((pp {K) +z p (pp{cr m + 5)) = {h n 0 p mp c m ) + (£(h n , c m ) ® p 5). 

From this we obtain 

pp~ rn 

p m / t + — (l — {h n 0 P {c m + 5)) 2 ) dh 
Jo ^ 

rp~ m 

1 + ~ (l — {{h n ©p mp C m ) + ( £,{h n , Cm) ©p h)) 2 ) df>. 


= p 


We now use j 6 j Lemma 3], which states that for any / G L 2 ([0,1]) and any y G [0,1), we 
have 


f{x) dx = f{x © p y) dx. 


(4) 


This yields 


P" 


1 + y (l — {h n 0 p {c m + <5)) 2 ) d<5 — 


= P" 


i o 

7i 


l + ^(l-((h n ©7 p a m ) + <f) 2 ) dS 


— 1 + y (1 - (M 

Furthermore, in a similar fashion, 


Nsmp 


CTr 


'1 2'j_ (h m sm P rr 1_ 

1 ' p m 2 n m p 2m 6 ' 


P“ 


1 + 7 i (1 — {{h n ©p {c m + 5))) d<5 — 
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— p m f 1 + 71 (l — {{h n ©p mp c m ) + {£{h n , c m ) © p <5))) dh 
Jo 

= P m [ P 1 + 71 (1 - {{K ©r p ©«) + 5)) 

Jo 

Finally, let us deal with the expression 
rp~ rn 

p m 1 + 71 min {1 - (h n ®p (cr m + h)), 1 - {h k ©p (cr m + 5))} d<5 (5) 

Jo 

with k ^ n. Note that, ask ^ n, we cannot have h n (B p {c m + 5) = h k ® p {a m + 5). Suppose 
that 

hn ©p (©n © ^) ^ h k ©p {c m © 5). (6) 

Using the notation introduced above, we can rewrite (J6J) as 

{h n ffip mp c m ) + {£{h n , c m ) ffip 5) < {h k ffip mp c m ) + {£{h k , c m ) ©p 5). 

Again, since k ^ n, we cannot have 

{hn ©r O = {hk ©r p o, 

as this would also imply £(/i n , cr m ) = £(/ifc, cr m ), and would so yield a contradiction to (ED- 
Furthermore, it cannot be the case that 


{hn ©r O > (/l fe ©7 P O, 

since £{h n ,c m ),f;{hk,c m ) G [0,p _m ), and so we would also end up with a contradiction 
to (ED- Therefore, we see that (ED automatically implies 

{K®^ Cm) < {hk®™* Cm)- ( 7 ) 

Suppose now, on the other hand, that (ED holds. Then, since £,{h n ,c m ),t;{h k ,c m ) G 
[0 ,p~ m ), also (J 6 D must hold. We have thus shown that (ED and (J2J) are equivalent. 
Suppose now in the analysis of (ED that (ED holds, i.e., 

rv~ m 

p' n / 1 + 71 min I 1 - {hn ©p {c m + 5)),1- {hk ©p {c m + 5))} d<5 = 

Jo 

rp~ m 

= p m 1 + 7i (1 {hk ©p {c m + $))) d<5 
Jo 

= P m [ P 1 + 71 (1 - {{hk ®r p O + m k , a m ) ffi p £))) d5. 

Jo 

Using the equivalence between (ED and ©, and again (ED, we see that the latter expression 
equals 


V 


m 



1 + 7! (min (1 - {h n ©p mp c m ), 1 


{h k ® s p mp c m )}-m k ,a m )® p 6)) dS 


7 


= p" 


1 + 7i (min {l - (/i n ©p mp 0 > 1 “ ( h k 


Nsmp 




»)}-£) CM 


7i 1 


= ~ — + l + 7i min {l - (h n ffip” p cr m ), 1 - (h k ®^ mp cr m )} . 

A similar argument holds if the converse of (J6J) holds. 

Putting all of these observations together, we obtain 

P rn / e AT,l(-^p,JV ©p (<?m + S),K l |7 i) d<5 = ^1 + 


N—l 


n =0 


7i 


~n E ( 1 +t i 1 - (ft ” ®r p ®») 2 ) - 




2 


p 2m 6 


1 N ~ 1 / 1 

^E(-lV +1+7I ~ 7lP ‘"® rP ' 7m) 


n=0 
N -1 


1 /I \ 

+VJ E -f = + 1 + 71 min {1 - (K ®r p *»), 1 - (ft* ®r p *»)} 

oo b — C\ \ r / 


n,k =0 
n^k 


> 


i + * 
3 


v-i 


iV 


n=0 

N-l 


E(' + t(i-(^ 


Nsmp 


>) 2 )-|-+;2(ft„®7-,7 m ) 71 1 


iV 2 

1 

'iV 2 


n=0 
IV-1 


^E p+71 u-K.® r> 


2 2p 

1 


2 4p : 


,2m 


2 p" 


^ f 1 + 7i min jl - (h n ffip mp a m + 

, \ v \ 


= i + 


n,k =0 
n^k 


7i 


N-l 


2 p 1 


1 A , 1 - f h k ffi! mp a m + 1 


2p r ‘ 


N-l 


N 


n =0 


E(i+|-( i -®®r 1 ®7 


Jp E (I + 71 min { I - Hi, 


l ) •> 1 (^/c 


^m)}) 


n,/c=0 


= e 2 N>1 (H PiN ffip a m , Ad, 71 ). 
The result follows. 


□ 


For two point sets X = {x 0 , xi, ..., aiAr_i} in [0,1 ) S1 and Y = {y 0 , y 1 ,.. ., 2 /jv_;l} 
in [0,1 ) S2 we write (A", Y ) to denote the point set consisting of the concatenated points 
(x k ,y k ) = Ofc,i,.. .,x kiS1 ,y ki i,... ,y k , S2 ) for k = 0,1,..., N - 1. 

Lemma 2. Let P S: n be a point set of N points in [0, l) s . Let H P} n be as in Lemma 0 and 
let m € N be minimal such that N < p m . Furthermore, let a m € Q(p m ). Then it is true 
that 


e 2 N 9 s+i((Ps,N, d c m ), K s+ln ) < p m 


^7\r,s+1 ((-^s ,Ni Hp N ffip (c m H - (5)), -K's+lj'y) ^5. 








Proof. The proof is similar to that of Lemma [Q 


□ 


3 The CBC construction 


In this section, we analyze the following CBC construction of a mid-simplified p-adic shift 
to obtain p-adically shifted Halton sequences with a low integration error. 

Throughout this section, let s, N G N be given and let p = (pi,---,p s ) £ with 
pairwise distinct components pj. For j £ [s] let m,j 6 Nbe minimal such that N < p'C. 
Let H p n be the point set consisting of the first N elements of H p . To stress the dependence 
of the worst-case error on the p-adic shift we write in the following 


e N ,s((r) ■= e N:S (H PyN ®p ld <r, K sn ) 


for cr G Qfp™ 1 ) X • • • X Q(p™ s ). 

We propose the following algorithm. 


Algorithm 1. (1) Choose a i G Q (p™ 1 ) to minimize e 2 N1 (a) as a function of a. 

(2) For 1 < d < s — 1, assume that oq,..., cr rf have already been found. Choose cr d +1 £ 

to minimize 

( 8 ) 

as a function of a. 

(3) If d < s — 1 increase d by 1 and go to Step 2, otherwise stop. 


Remark 1. We remark that Algorithm [T] makes the main result in [6] much more explicit, 
as the algorithm only needs to check a countable number of possible candidates for the 
p-adic shift. A slight drawback of our method is that the effective CBC construction of 
good p-shifts has a cost of 0(s 2 N 3 ) operations, which is still large. Further improvements 
with respect to the construction cost are a demanding problem for future research. 


The following theorem states that Algorithm [Tj yields p-adically shifted Halton se¬ 
quences with a low integration error. Note that the error bound is of the same order as 
the one in Theorem [EJ 

Theorem 2. Let the notation be as above, and let d G [s]. Assume that cr s = (ay,..., a s ) 
has been constructed according to Algorithm. QJ Let <Td ■= ( 01 ,, o~d)- Then 

e 2 N4 {(T d ) < ^ (n f 1 + 2 7j( io g^)^-) + n( i +7,)n( i +^)) • ^ 

\j=1 v 613 ' j= 1 3=1 / 

Proof. We show the result by induction on d. For d — 1 we have 

e N AH pl , N ffipi 0 ) A 1 . 71 ) dcr 

H P i,N ©pi d(f 


1 P LW r& 

— E p™‘ / 

Pi Jl/ri 


(t+A/rT 1 


P 2 

e NA 
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where we applied Lemma [D Hence there exists a a[ E QQo™ 1 ) such that 


e 2 N, i( a i) < / e N,i( H Pi,N ® P1 <J,K lni )da 


< i J (l + 2 7 l (logW) I ^_)+(l +7l )(l + ^fl), 


where we used m Theorem 1] for the second inequality. Since op is chosen by Algorithm Q] 
to minimize e 2 N1 (a), it follows that the result holds for d — 1 . 

Suppose the result has already been shown for some fixed d E [s — 1]. Assume that 
<t ( i — (< 7 i,..., ad) has been obtained by the CBC algorithm. Since a^+i is chosen in order 
to minimize the squared error (l 8 j) . we have (where we write with some abuse of notation 
(cr<2, &d+ 1) := (o"i, • • •, &di &d+ 1)) 


1 


Pd+i 


e N,d+l(( Cr d, ©Z+l)) A m d+1 

Pd+1 v =0 


y e N,d+l ( ( a d, 


™d +1 

Pd+1 


Using Lemma El we now see that, for any v E {0,... EPd+V ~ 1}, 


v 

e ~N,d+l 1 1 a 'di m d+ 1 

Pd+1 

— Pd+1 j e N,d +1 ^ (^Hp d ,N ©p <^di Hp d+1 ,N ©p d+ i ^m d+1 + ^ ) ) j Kd+ 1,7 ) d(5, 


where p d := {p \.... ,Pd), and hence 


e N,d+i((<7d , (?d+ 1 )) < / e N,d+i(( H Pd,N ©p R tr d , Hp d+ 1 ,JV ©Pd+1 v), K d+ 1 , 7 ) da. 

Jo 

We denote the points of H Pd ,N®f d (Td by x n (^n,i; ■ ■ ■ 5 %n,d)i and the points of y 

by h n . Due to (J2]) , we obtain 


d+1 


N -1 


j e N,d+i((H Pd} N ©p ^mi Hp d+1 ,N ©p d+1 a), Arf+i j7 ) da 

n i 1 + ^C 1 ~ X nj)) / ( 1 + ^y 1 ( 1 -( /l n©p d+1 ^) 2 )) da 

J=1 J 

' d 

J (1 + lj min{l - x nJ , 1 - x kJ }) 


- V 

AT 

n —0 
1V-1 

iW S 

n,/c=0 Lj=l 


X / (l + 7d+i min{l - (h n © Pd+1 


a),l - (h fc ffip d+1 a)}) da. 
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Let now 


h ■= 



1 + ~ V ln ®Pd+i a ) 2 )) dcr, 


and 


h = (l + 7d+i minjl - (h r 


7 Pd +1 


a),l - (h k ® Pd+1 a)}) da. 


Using (gD, we obtain 

h = 

Let us now deal with J 2 



ld+1 


1 + ^(1-a 2 ) da = 1 + 


7d+l 


e=o 


Pd+l,ld+l 


(£)(3e(h n )/3e{h k ), 


where for l = £ a -iPd+i + • • • + £iPd+i + £o with £ a _i ^ 0 we have 


1 + 

' Pd+l,7d+l — 'l ld+1 ( 1 


2 Pd+l \ sin 2 a _ i 7T /pd+ 1 ) 3 


if t = 0, 
if £ ± 0. 


Altogether, we obtain 
e Ar,d+l(( cr rfi °d+l)) 


< 2+1 


< 


1=1 


nb+i)-vE 


N—l 


n =0 Li=l 




X. 


n,]j 


1 + 


ld+1 


N 2 


N-l 

E 

n,k =0 b'=l 


J (1 + 7i min{l - x nJ , 1 - x kJ }) 


E 

£=0 


Pd+li7d+l 


(£)l3e(h n )f3i(h k ) 


= 1 + 


ld+1 


N-l d 


n(^i)-^Eno + |d- 


X 


n,3 


J = 1 


n=0 J 1 


1V-1 d 


Jp Y II + 7 * min {! - ^nj, 1 - 


n,k=0 j=1 


N-l / <2 


IV 2 
1 + 


Y ( n ^ + U mini 1 - Znj, 1 - Xfcj}) ) Y 

7=1 

&N,d{ (T d) + U, 


Pd+l>7d+l 


(£)/3e(h n )/3 e {h k ) 


n,k= 0 \j=l 

7d+l \ 2 


( 10 ) 


where 


JV-1 / d 


T:= J^Y1 ( minjl - x nJ , 1 - x Kj }) ) Y 


Pd+i,id+i 


(£)/3i(h n )/3 e (h k ). 


n,k =0 \j = l 


1=1 
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Since min{l — x n j, 1 — Xkj} < 1 we obviously have 


r< (nii + ^lE 

\j=1 


Pd+l^d+l 




f=\ 


N—l 


Jf £ Mhn) 


n —0 


(ii) 


From the proof of P Theorem 1], it can easily be derived that 

2 


£ 

e =i 


Pd+i,id+i ' 


N—l 


N 


£ A (hr. 


n =0 


< 1 id+igPd+i + 7d+i £ + jd+i 


N 2 


6 P 9 d+ i 


for arbitrarily chosen g G N. By choosing g = |_21og pd+i IV J and inserting into (ITT]) . we 


arrive at 


T < 


< 


N 2 

1 

N 2 


3 = 1 


1(1 + 7i) ( ( 7d+i (log AT) 




Pd+i 

logPd+i 


+ 


Id+iPd+i 


6 


1 + 


7d+i 


(7d+i(log N) 




Pd+i 

logPd+i 


Jd+iPd+i 

6 


d +1 


i=i 

d 


1(1 + 2 Tj (logJV) 


Pi 


log Pj 


IF+a)!! 1 + 


7 j n 

6 


i=i j= i 

On the other hand, we have, using the induction assumption, 
7d+l \ 2 


1 + 


< 1 


< 


e N,d( a d) 

7d+l \ J_ 
N 2 


vi=i 


J(l + 2 7 j (logiV) 


Pi 


log p. 




i 

N 2 


1 + 7 J+i(logJV) kStr)n 


3 S j =1 j =1 

d / 2 

l + 2 7i (log IV) 


7iPi 


Pj 


log Pj 


d+l 


id+aoii o+ 


7 jPj 
6 


3=1 3=1 

Combining equations (TT 2 |) and (fT3j) . and inserting into (fTUl) . we obtain 


1 


f d+l 


■'iVjd+l' 


((o-d^d+O) < ^2 ( II ( 1 + 2 7j(logT) 




log Pj 


0=1 X 3=1 

This is the result for d+l, and the theorem is shown. 


d -\-1 d -\-1 

II( 1 + a)]I ( 1 + 

3 = 1 


7 jPj 


( 12 ) 


(13) 


□ 


References 

[1] J. Dick, F. Y. Kuo, F. Pillichshammer, I. H. Sloan: Construction algorithms for 
polynomial lattice rules for multivariate integration, Math. Comp. 74, 1895-1921, 
2005. 


12 

























[2] J. Dick, F. Pillichshammer: Multivariate integration in weighted Hilbert spaces based 
on Walsh functions and weighted Sobolev spaces, J. Complexity 21, 149-195, 2005. 

[3] J. Dick, F. Pillichshammer: Digital Nets and Sequences. Discrepancy Theory and 
Quasi-Monte Carlo Integration , Cambridge University Press, 2010. 

[4] M. Drmota, R. F. Tichy: Sequences, Discrepancies and Applications , Springer, 
Berlin, 1997. 

[5] J.H. Halton: On the efficiency of certain quasi-random sequences of points in evalu¬ 
ating multi-dimensional integrals. Numer. Math. 2, 84-90, 1960. 

[6] P. Hellekalek, P. Kritzer, F. Pillichshammer: Open type quasi-Monte Carlo integra¬ 
tion based on Halton sequences in weighted Sobolev spaces. Submitted, 2014. 

[7] Korobov, N.M.: Number-theoretic methods in approximate analysis. Gosudarstv. Iz- 
dat. Fiz.-Mat. Lit., Moscow, 1963. (in Russian) 

[8] P. Kritzer, F. Pillichshammer: On the component by component construction of 
polynomial lattice point sets for numerical integration in weighted Sobolev spaces. 
Unif. Distrib. Theory 6, 79-100, 2011. 

[9] L. Kuipers, H. Niederreiter: Uniform Distribution of Sequences, John Wiley, New 
York, 1974. 

[10] F. Y. Kuo: Component-by-component constructions achieve the optimal rate of con¬ 
vergence for multivariate integration in weighted Korobov and Sobolev spaces, J. 
Complexity 19, 301-320, 2003. 

[11] H. Niederreiter: Random Number Generation and Quasi-Monte Carlo Methods. No. 
63 in CBMS-NSF Series in Applied Mathematics. SIAM, Philadelphia, 1992. 

[12] E. Novak, H. Wozniakowski: Tractability of Multivariate Problems. Volume I: Linear 
Information, EMS Tracts in Mathematics, 6. European Mathematical Society (EMS), 
Zurich, 2008. 

[13] I.H. Sloan, S. Joe: Lattice Methods for Multiple Integration. Oxford University Press, 
New York and Oxford, 1994. 

[14] Sloan, I.H., and Reztsov, A.V.: Component-by-component construction of good lat¬ 
tice rules. Math. Comp. 71, 263-273, 2002. 

[15] I.H. Sloan, H. Wozniakowski: When are quasi-Monte Carlo algorithms efficient for 
high-dimensional integrals? J. Complexity 14, 1-33, 1998. 

[16] X. Wang. A constructive approach to strong tractability using quasi-Monte Carlo 
algorithms. J. Complexity 18, 683-701, 2002. 

Authors’ addresses: Peter Kritzer, Friedrich Pillichshammer 

Department of Financial Mathematics and Applied Number Theory 

Johannes Kepler University Linz 

Altenbergerstr. 69, 4040 Linz, Austria 

E-mail: peter.kritzer®jku.at, friedrich.pillichshammer@jku.at. 


13 


